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ABSTRACT 

We present a Bayesian reconstruction method which maps a galaxy distribution from 
redshift-space to real-space inferring the distances of the individual galaxies. The 
method is based on sampling density fields assuming a lognormal prior with a likeli¬ 
hood given by the negative binomial distribution function modelling stochastic bias. 
We assume a deterministic bias given by a power law relating the dark matter density 
field to the expected halo or galaxy field. Coherent redshift-space distortions are cor¬ 
rected in a Gibbs-sampling procedure by moving the galaxies from redshift-space to 
real-space according to the peculiar motions derived from the recovered density field 
using linear theory with the option to include tidal field corrections from second order 
Lagrangian perturbation theory. The virialised distortions are corrected by sampling 
candidate real-space positions (being in the neighbourhood of the observations along 
the line of sight), which are compatible with the bulk flow corrected redshift-space 
position adding a random dispersion term in high density collapsed regions. The lat¬ 
ter are defined according to the eigenvalues of the Hessian. This approach presents an 
alternative method to estimate the distances to galaxies using the three dimensional 
spatial information, and assuming isotropy. Hence the number of applications is very 
broad. In this work we show the potential of this method to constrain the growth rate 
up to k ~ 0.3 h Mpc -1 . Furthermore it could be useful to correct for photo-metric 
redshift errors, and to obtain improved BAO reconstructions. 

Key words: cosmology: large-scale structure of the Universe - cosmology: theory - 
galaxies: general - methods: observational - methods: numerical 


1 INTRODUCTION 

Galaxy redshift surveys produce the three-dimensional dis¬ 
tribution of luminous sources tracing the underlying dark 
matter field. However, their infered line-of-sight position is 
a combination the so-called Hubble flow, i.e. their real dis¬ 
tance, and their peculiar motion. The modifications pro¬ 
duced by this are referred to as redshift space distortions 
(RSD). Many astronomical studies are limited by 
distortions, such as a proper environmental study 


these 

(Nuza 
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|et al.||20l4 ). Nevertheless, RSD can also be used to con¬ 
strain the nature of gravity and cosmological parameters 
(see e.g. Berlind, Narayanan &; Weinberg] |2001 | Zhang et 


aL][2007| Jain &; Zhang||2008| |Guzzo et al.|| 2008 Neseris^ 

Perivolaropoulos|2008| |Song &; Koyama|2009| |Song &; Perci- 

val||2009| |Percival &; White||2008l |McDonald Seljak||2009| 

White, Song &; Percival||2009| |Song et al.||2010| |Zhao et al.| 


2010| Song et al.|2011[ for recent studies). The measurement 
of RSD have in fact become a common technique (Cole, 


Fisher & Weinberg 1995; Peacock et al.|2001; Percival et al. 

2004j da Angela et al. 

2008; Okumura et al. ;2008J^ Guzzo 

et al. 2008; Blake et a 

. 2011; Jennings, Baugh & Pascoli 
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Raccanelli|2012| |Reid et al.|2012| |Okumura, Seljak &; Desc- 

jacques|2012||Samushia et al.|2013||Zheng et al.|2013| Blake 


e^aLj 2013JJdel^Torr^et^LJj2013J^ Samushia^^il 


Sanchez et al.| 

2014 

Okumura et al. 

|2014 


2014 


usually based on the large-scale anisotropic clustering dis¬ 
played by the galaxy distribution in redshift-space, although 
iV-body based models for fitting the data to smaller scales 


have been presented in Reid et al. (2014). 


In this work we use the full three dimensional distri¬ 
bution of galaxies and correct in a statistical Bayesian way 
their individual positions according to a physical model de¬ 
scribing the relation between the density and peculiar veloc¬ 
ity fields which depends on the growth rate. This relation 
becomes more complex if one wants to correct for dispersed 
RSD as we will discuss below. The main finding of this work 
is that already very simple assumptions, such as a lognormal 
prior for the density field, a power-law galaxy bias, and lin¬ 
ear theory to infer the peculiar motions with an adjustable 
smoothing scale, can yield accurate reconstructions of the 
real-space position of galaxies beyond the Kaiser-factor. 

In the following section we will present our methodology 
and subsequently we will present the application on a mock 
galaxy catalogue. Finally we will summarise our findings and 
present our conclusions. 


2 METHOD 


Our method essentially follows the algorithm (ARGO-code) 
proposed in Kitaura & EnBlin (2008) and Kitaura, Gall erani| 


&; Ferrara (2012) by iteratively sampling the density and 


peculiar velocity fields within a Gibbs sampling process (for 
the pioneering iterative correction of linear RSD see |Yahil| 
et al.j 1991): 


S ^ V(6\N({r}),{b p },C,w) , (1) 

M ^ v ({r}\{s ohs },{v(8,M,r s ,-H,{v P })}) ■ (2) 

In the first Gibbs-sampling step: 5 is the dark mat¬ 
ter over-density field, V (5 \ IV({r}), { 6 P } , C, w) the poste¬ 
rior distribution function of density fields given the num¬ 
ber counts of galaxies in real-space on a grid AT({r}), a co- 
variance matrix C = (S5^) (the power-spectrum in Fourier- 
space), a set of parameters describing galaxy bias {5 P }, and 
a three dimensional completeness w. In the second Gibbs- 
sampling step we obtain the real-space position {r} from the 
redshift-space position of galaxies and the sampled peculiar 
velocity fields {v}, which depend on the large scales, on the 
over-density field, the growth factor /q, and a smoothing 
scale rs; and in the nonlinear regime additional information 
about the Hessian PL and parameters describing the velocity 
of galaxies {u p } . 


2.1 Density field reconstruction 

We will restrict ourselves to a simple set of assumptions 
and demonstrate that they are sufficient for the scope set in 
this work, namely correct for redshift-space distortions. The 
first Gibbs-sampling iteration assumes that the galaxies are 



Figure 1. Upper panel: power spectra of the catalog in redshift- 
space (red), in real-space (black); and of the reconstructed cata¬ 
logues with only coherent flows (dashed green) with 1-sigma con¬ 
tour based on 1000 reconstructions (magenta), coherent flows in¬ 
cluding virialisation corrections (dashed blue) with 1-sigma con¬ 
tour (cyan). Lower panel: ratio between the power spectra in 
redshift-space P s (k) and in real-space P r (k) with the same colour 
code. 


in real-space and applies a Hamiltonian sampling technique 
( Jasche k, Kitaura||2010 ), to obtain the density field within 


the Lognormal approximation with stochastic bias (Kitaura 
et ah]2010). In particular we consider the negative binomial 


distribution function with nonlinear bias, as introduced in 


At a, Kitaura &; Muller (2015). We note that the prior could 


be improved introducing higher order correlation functions 
(Kitaura ||2012 ) or perturbation theory based methods (see 


e.g. Kitaura 2013). Nevertheless, we want to remain with 
simple models which require the last number of assumptions 
and parameters. 

The deterministic nonlinear bias we consider here is re¬ 
stricted to the following form: 


Pg — 7Pm , 


(3) 


with ph being the expected number density of galaxies in 
a cell within a volume embedded in a mesh, pm being the 
dark matter density on the same mesh, a being the power 
law bias, and 7 the proportionality factor dependent on 
the mean number density. Here we have neglected threshold 
bias, which would yield a more precise description in terms 
of higher order statistics ( Kitaura et al.|20l4 Kitaura et al. 
2014b). We do not consider non-local bias, which will be 


introduced in a forthcoming publication. 


2.2 Peculiar motion reconstruction 

The redshift-space position of a galaxy s obs is composed 
by its real-space position r and the coherent and dispersed 
components t;£ oh (r) and v%(r), respectively: 

s ohs = r + Vr° h (r) + Vr(r ), (4) 

with the subscript r denoting the projection along the line- 
of-sight. Let us discuss both redshift-space contributions 
separately below. 
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Figure 2. Left panel: Correlation function of the catalog in redshift-space £ s (r) (red), in real-space £ r (r) (black); and of the reconstructed 
catalogues with coherent flows including virialisation corrections (dashed blue) with 1-sigma contour of 20 samples (cyan). Additionally, 
the reconstructed catalogues with coherent flows corrections only are shown: magenta without Gaussian smoothing, green with optimal 
Gaussian smoothing of rs =7 h -1 Mpc. Right panel: corresponding quadrupoles Q s (r ), Q r (r). In the lower panels the corresponding 
ratios with respect to the signal in real-space are shown. 


2.2.1 Coherent redshift-space distortions 

Coherent redshift-space distortions are responsible for the 
squashing effect of galaxy clusters along the line of sight (see 
Hamilton 1998 for a review). This produces an enhancement 


of power on large scales, the so-called Kaiser factor (Kaiser 


1987). 


For simplicity we focus in this work on linear theory 
based on the density field (J). We are aware that improve¬ 
ments could be found based on the linearised density field 
(log(l + (5)), which we get for free in our formalism. How¬ 
ever, it is also true that the logarithmic transformation in¬ 
troduces a constant offset which we do not want to consider 
here, as it was shown in |Neyrinck et ah (2009) (for more 
general relations including linearisations see | Kit aura et ah] 
2012 and references therein). We simply consider a Gaus¬ 


sian smoothing (with radius rs) of the overdensity field to 
optimise the velocity divergence to overdensity relation. For 
this component we do not assume any stochastic component, 
but directly move each galaxy from its redshift-space posi¬ 
tion s obs to its coherent real-space position r z coh at iteration 
i, according to the coherent bulk flow motion Vr° h,z (r' l ~ 1 ) 
based on the real-space position from the previous iteration 
z — 1: 


i obs „ coh,z ( i— 1\ 

^coh ® V r V ) ; 


(5) 


with v r = (v-r)r/(Ha), where v is the full three dimensional 
velocity field, r is the unit sight line vector, H the Hubble 
constant and a the scale factor. 


2.2.2 Dispersed redshift-space distortions 

Further in the nonlinear regime in deep gravitational wells, 
galaxy clusters are quasi-virialised, and produced so-called 
fingers-of-god which are elongated structures along the line- 
of-sight (Jackson 1972). This reduces the clustering towards 
small scales. A method aiming at correcting these effects 
needs thus to enhance the clustering of galaxies in clus¬ 
ters. We do so by first making sure that the galaxies come 


from collapsed regions in real-space. To this end we demand 
that the real-space candidate positions along the line-of- 
sight come from cells in which all the eigenvalues of the 
Hessian are positive (Hahn et al. 2007) and above a certain 
overdensity threshold J t h- This supposes another quantita¬ 
tive application of the cosmic web classification, in addition 


to the first of such kind presented in Zhao et al. (2015). For 
each real-space candidate position r k we sample the viri¬ 
alised motion component from a Gaussian with dispersion 
= rjPmirk) (r): 

vf(r k ) = Q (vf(r k ) | cr(r k ), \/X(H(5) > 0)) r , (6) 


with additional parameters p and e (see |Kitaura et al.' 2014|). 

We can then find coherent real-space positions r kco h = 
r k + vf(r k ) and select the closest one in each iteration: 
min(|r co h — T*fccoh|)- Additionally, to ensure that the cluster¬ 
ing of clusters is enhanced we demand that the local density 
of the candidate’s real-space position is larger or equal than 
in the previous iteration: 5(r k ) x ^ J(r/e)* _1 . We note that a 


similar approach was taken using the KlGEN-code (Kitaura 
2013) to correct for fingers-of-god in a forward approach (see 
Hefi, Kitaura fe Gottlober| 2013). 


3 APPLICATION TO MOCK DATA 

We will focus in this work on mock galaxies constructed 
using the halo abundance matching technique based on the 
CMASS LRG sample (see Rodrfguez-Torres et al in prep), 
and the BigMultiDark (BigMD) iV-body simulation. 


3.1 Reference iV-body simulation and CMASS 
LRG galaxy catalogue 

In particular, we performed a halo abundance matching 
using the redshift 0.5763 from one of the BigMD simula- 
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Figure 3. The upper panels show the 2D power spectra P(k±,k\\) corresponding to the mock galaxy catalog in real-space (colour- 
coded contour regions), in redshift-space (dashed lines), including the reconstructed galaxy field in real-space (solid lines): without any 
additional smoothing (2nd panel), with an optimal smoothing of rg = 7h -1 Mpc (3rd panel), including virialised RSD correction (4th 
panel). The lower panels show the corresponding ratio of the 2D power spectra, (left panel:) between the true catalogue in redshift-space 
and in real-space, and (for the 3 panels on the right) between the corresponding reconstructed real-space catalogue and the true catalogue 
in real-space. 


tion^ ( Klypin et al.||2014 ), whic h was performe d using the 
TreePM N -body code gadget-2 ( Springel|2005 ) with 3840 s 
particles and the volume of (2.5 h~ 1 G pc) , in a framework 
of Planck ACDM cosmology with {Q m = 0.307115, Qb = 
0.048206, ag = 0.8288, n s — 0.96}, and the Hubble param¬ 
eter (Ho = 100 /ikms _1 Mpc _1 ) given by h — 0.6777. We 
used the spherical overdensity (BDM) halo catalog. The first 
step to generate the galaxy catalog is to modify the maxi¬ 
mum circular velocity (Wax) of each objects adding a gaus- 


sian noise: 


^ e a r-Wax(l+AT(0,a)) (7) 

where Af(0, cr) is a gaussian random number with mean 0, 
and standard deviation a. Then, we sorted all objects by 
Wax and then, we selected objects starting from the one 
with larger Wax and we continue until we get a number 
density of 3.29x 10 _4 /i 3 /Mpc 3 . Finally, we fit the clustering 
of the BOSS CMASS sample using the scatter parameter a. 
We consider for our study the catalogue in a sub-volume of 
(1.25 /i -1 Gpc) 3 . 


3.2 Reconstruction results 

To avoid systematic effects in our study, we neglect incom¬ 
pleteness due to the survey geometry or the radial selec¬ 
tion function, and neglect light-cone evolution effects. All 
these issues will be considered in a forthcoming work (Ata 


1 http://www.multidark.org/MultiDark/ 


et al in prep). We compute the redshift-space position for 
each galaxy in the plane parallel approximation. We then 
compute the corresponding density field on a cubical region 
with 128 3 cells. This permits us to make efficient compu¬ 
tations of the 2D power spectrum. We note however, that 
our ARGO code has been implemented to deal with non¬ 
parallel RSD. We obtain converged Gibbs-sampling chains 
after about 1000 iterations in terms of converged power spec¬ 
tra and also according to our convergence study demon¬ 
strated in (Ata, Kitaura & Muller 2015). We have tested 
here different levels of deviation from Poissonity. As in |Ata,| 


Kitaura & Muller (2015) we find that the stochastic bias 


can enhance the power towards high fcs, however, this is 
not essential to our work. The smoothing scale rs and the 
threshold Jth can be tuned to compensate for that. Since 
we do not aim at getting the perfectly unbiased dark mat¬ 
ter density field in this work, but to correct for RSD, we 
will show only results with very high /3, effectively sampling 
from the Poisson likelihood. An extension of this work in¬ 
vestigating also the reconstructed dark matter field will be 
presented in Ata et al. in prep. 

We run three reconstruction chains, two correcting only 
for coherent RSDs, and the third one correcting also (par¬ 
tially) for virialised RSDs. 


(i) The first RSD correction is based on coherent peculiar 
motions directly derived from the density field on a mesh 
using linear theory. 

(ii) The second one uses a smooth density field with Gaus¬ 
sian smoothing radius of 7 h~ x Mpc obtained in a parameter 
study to optimally correct for RSD up to higher ks. 
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(iii) The third one includes virialised corrections as de¬ 
scribed in §2.2.2| 

The resulting power spectra considering 1000 recon¬ 
structions only for the latter two cases for clarity are shown 
in Fig. [I] The first case yields a similar result on large scales, 
however underestimating the power in intermediate and 
small scales. Here we can clearly see that the Kaiser factor is 
corrected and that the nonlinear RSD correction increases 
the clustering power towards small scales being closer to 
the true real-space catalogue. Fig. [2] shows the monopole 
and quadrupole in configuration space. We can see from 
these plots that the two different reconstructions: coher¬ 
ent with optimal smoothing or coherent and virialised mo¬ 
tions corrections yield similar results. The left panel shows 
how accurately the real-space BAO can be obtained from 
redshift-space. However the right choice of the smoothing 
scale is important here, or a nonlinear virialised treatment 
to obtain a close BAO peak to the true one. We will quan¬ 
tify this improvement investigating also BAO reconstruc¬ 
tion in a forthcoming work. Interestingly, the full nonlinear 
RSD correction algorithm shows better agreement with the 
true real-space quadrupole not only on small scales being 
closer to zero, but also at large scales, displaying less arti¬ 
ficial spikes, present in the pure coherent RSD corrections. 
However, the different quality in the reconstructions can be 
better appreciated in Fig. [3] showing the 2D power spec¬ 
tra. The anisotropic pattern in redshift-space can be clearly 
seen. In particular, the enhancement of power due to the 
Kaiser effect is very prominent. A reconstruction of the pe¬ 
culiar velocity field with our method corrects the RSD in 
the catalogue in a remarkable way. We can see that an opti¬ 
mal choice of the smoothing scale can considerably improve 
the reconstruction. This is furthermore improved to scales 
of about k ~ 0.3 h Mpc -1 when including virialised motions 
corrections. 


4 SUMMARY AND DISCUSSION 


We have presented in this work a Bayesian technique to 
correct for both coherent and virialised redshift-space dis¬ 
tortions present in galaxy catalogues by estimating the dis¬ 
tance to the individual galaxies. We have demonstrated that 
this technique is accurate at least up to k ~ 0.3 h Mpc -1 
in the isotropisation of the 2D power spectrum based on 
precise galaxy mock catalogues describing the CMASS LRG 
sample. Although the method is general enough there to be 
precise down to far smaller scales, as indicated by the recov¬ 
ered power spectra. 

An application of this technique to the BOSS DR12 
data including power spectrum sampling will be presented 
in a subsequent publication (Ata et al. in prep). 

While traditional RSD measurements focus on the 
growth rate, an approach like the one presented in this work 
is complementary and more general. We refer to a recent ap¬ 
plication of a joint analysis of density fields and anisotropic 
power spectra including growth rate estimation, see Gr anett| 
et al. (2015). The advantage of the approach presented in 


the present work is that it deals with nonlinear structure for¬ 
mation, nonlinear and stochastic galaxy bias, yielding also, 
as a by-product the real-space positions of the individual 
galaxies. 


This technique is promising for a broad number of ap¬ 
plications, such as correcting for photo-metric redshift-space 
distortions including the cosmic web information, or to make 
precise environmental studies, as demonstrated in N uza et| 


al. (2014) with a similar forward method recovering the cor¬ 


responding primordial fluctuations. We have demonstrated 
in particular that it is a potentially interesting technique for 
the estimation of the growth rate, or for an improved BAO 
reconstruction. These topics will be investigated in detail in 
forthcoming publications using this technique. 
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